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Abstract 

We study the vortex dynamics in an evolutive flow. We carry out 
the statistical analysis of the resulting time series by means of the joint 
use of a compression and an entropy diffusion method. This approach 
to complexity makes it possible for us to establish that the time series 
emerging from the wavelet analysis of the vortex dynamics is a weakly 
chaotic process driven by a hidden dynamic process. The complexity 
of the hidden driving process is shown to be located at the border 
between the stationary and non-stationary state. 

1 Introduction 

In the last few years an increasing interest has been devoted to the study 
of weak chaos, which has become for many authors the paradigm itself of 



1 



complexity. The characterization of weak chaos, as it emerges in the context 
of deterministic chaos from the work of Zaslavsky et ah Ij, imphes the 
existence of a phase space where regular regions coexist with random regions. 
Some authors, see for instance Ref. 0, propose a dynamic condition where 
randomness and order are simultaneously present. The latter is a property of 
systems, usually organized systems, where the clock and coin flip paradigms 
apply at the same time [S]. Indeed, it is plausible that a condition of total 
randomness is incompatible with the self-organization that is expected to 
play a key role in such dynamical systems jl]. 

In this sense, dynamic models of the same kind as the Manneville map jS] 
might be considered to be representative of the complexity of living systems, 
as suggested, for instance, in the work of Refs. [HI Hj . The Manneville map 
has been discussed as a model of intermittency years ago by Gaspard and 
Wang [HI and more recently rigorous proofs and improvements of the above 
results have been obtained by Bonanno ^ and Galatolo JTOJ. 

We recall that the Manneville map f{x) = x+x^ (mod. 1) is characterized 
by a parameter z > 1, which controls the frequency of the random events. 
We can consider the waiting time between two consecutive random events as 
time of sojourn in a condition of order. The corresponding distribution of 
waiting times is an inverse power law with index n = z/{z — l). The condition 
2; = 2 is the border between the region where the waiting time distribution 
has a finite mean [z < 2) and the region where this mean waiting time is 
infinite {z > 2). This latter behavior has been named sporadic dynamics [S]. 

From a statistical point of view, the condition z < 2 does not mean ordi- 
nary statistical mechanics. In fact, the region | < 2; < 2 is characterized by 
waiting time distribution with divergent second moment and for this reason 
it was shown to be the dynamic source of Levy statistics ^] . This form 
of statistics led the authors of Ref. ^H] to coin the term of strange kinetics to 
point out the anomalous nature of the diffusion process emerging from this 
condition. The same authors considered the strange kinetics as a synonym 
of weak chaos. 

Moreover, for a better understanding of the balance between order and 
randomness in the Manneville dynamical systems, it is suitable to take 
also the Kolmogorov- Sinai (KS) entropy ^3] into account. In the region 
1 < z < 2, the KS entropy is positive. The entropy vanishes when the 
parameter z enters the region z > 2, therefore generating a non-stationary 
condition. However, the characterization of weak chaos by means of non- 
ordinary statistics is not sufficiently fine for our purposes, since it applies 
to both the case when the KS entropy is finite and to the case when the 
KS entropy vanishes. Hereafter, we shall adopt the following definition of 
weakly chaotic systems: systems whose KS entropy is zero but which are not 
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periodic. Therefore, following our definition, the Manneville map is weakly 
chaotic only when the parameter z lies in the region z > 2. 

Nevertheless, since our definition of weakly chaotic dynamical system is 
focused on a vanishing KS entropy, we have to use an even finer indicator to 
classify the different types of weak dynamics. For instance, the logistic map 
at the Feigenbaum point represents a weakly chaotic scenario that arises 
from a period-doubling cascade. The authors of [15, 16j showed that the 
logistic map at the chaos threshold turns out to be mildly chaotic and is 
not an intermittent form of weak chaos like the sporadic randomness defined 
by the Manneville-like systems with z > 2. In order to sketch the main 
differences between the two kinds of dynamics, we can consider the quantity 
of information which is necessary to describe n steps of a typical orbit. In the 
case of the intermittent Manneville-like systems with z > 2, the information 
grows as a power law n" (where the exponent a is 1/(2 — 1)), while in the 
logistic map at the chaos threshold it grows as the logarithm of n (see Section 
4andRefs. P [13 El CUl). 

The aim of this paper is to study a fluid dynamical process in the outer 
shear layer of a coaxial jet flow: to be precise, the roll-up and pairing of 
vortices. We reach the conclusion that the observed dynamics is an example 
of weak chaos, driven by a hidden intermittent process very well described by 
an idealized version of the Manneville model. In order to obtain our result we 
used three distinct techniques. The first one is the wavelet-Hilbert analysis 
[TS] . The second one is the Diffusion Entropy (DE) method [7] and, 
finally, the third one is the Computable Information Content (CIC) method 
of Refs. fK)\ IT7] . This is the second example of a joint application of the 
DE and CIC methods. The first is given by the work of Ref. (2j, where the 
reader can find a detailed discussion about the benefit stemming from the 
joint use of these two methods. 

As to the experimental apparatus used to acquire the data to study, the 
velocity signals are obtained in a coaxial jet flow (see the scheme in Figure^), 
described in more detail by Buresti et al. j21j. The analyzed configuration 
has an inner to outer jet diameter ratio Do/Di = 2.0, {Di = 50 mm), a 
velocity ratio Uo/Ui=1.5 with Uo = Q m/s and wall thickness L =5 mm. 
The axial and radial velocity components, u and f , were measured by means 
of X hot-wires placed inside the inner and the outer shear layers, at several 
axial positions downstream the jet outlet. The measurements at each point 
were carried out three times at a 5000 Hz frequency for 20 seconds. 

The outline of this paper is as follows. To make it as selfcontained as 
possible, in Sections 2, 3 and 4 we make a short review of the wavelet- 
Hilbert analysis, of the DE and of the CIC methods, respectively. In Section 
5 we illustrate the results of these techniques applied to the vortex dynamics 
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under study in this paper. Finally, in Section 6 we discuss the results that 
were obtained. 




Figure 1: illustration of the coaxial jet. 



2 The Wavelet-Hilbert analysis 

The experimental velocity signals acquired in many flow fields are characte- 
rized by fiuctuations induced by the passage or by the oscillation of vortical 
structures. For instance, in the developing shear layers of a coaxial jet the 
velocity signals are of increasing complexity with increasing distance from the 
jet outlet, and multiple dominating frequencies may be present, correspond- 
ing to fiuctuations with different degrees of modulation both in amplitude and 
in frequency. These fiuctuations are connected with the passage of vortical 
structures produced by the instability of the shear layers, and with their sub- 
sequent dynamics as well. Therefore, signal analysis procedures based only 
on conventional Fourier methods become largely inappropriate and differ- 
ent techniques, providing the time variation of the frequency and amplitude 
of the components present in the signals, should be used. A classical de- 
modulation technique, on which the analysis procedures used in the present 
investigation are founded, will now be briefly described. Any time-varying 
real signal, x{t) may be represented in the form 

x{t) = A{t) cos{(p{t)) . 
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This representation is not unique, but a so-called canonical pair may be de- 
fined for the quantities {A, ip) if we introduce the associated complex analytic 
signal Zx(t) defined by 



x{t) + iXHif) , 



where xuit) is the Hilbert tranform of x{t), given by 



XHit) = H[x{t)] = -P 




t-T 



dr . 



The canonical representation of x{t) is then 

x{t) = A^{t) cos[v9^(t)], 
hereby yielding the following definition of instantaneous frequency of x{t): 



This definition, which is perfectly correct from a mathematical point of view, 
acquires also a significant physical meaning when the analyzed signal is 
asymptotic, i.e., a sinusoidal signal which is slowly modulated in amplitude 
and frequency, x{t) = A(t) cos[u{t) ■ t + ipo]. In this case, provided the mod- 
ulation frequencies are sufficiently lower than the fundamental frequency of 
the signal, the time variation of amplitude and frequency may be directly 
recovered from the time variation of the modulus and of the phase derivative 
of the associated analytic signal. 

However, this procedure is applicable only if the signal is composed of a 
single asymptotic component; conversely, if multiple components are present, 
other techniques must be devised. One of these |22j is based on the appli- 
cation of the continuous wavelet transform and exploits the strict analogy 
existing between the wavelet transform of a signal and the extraction of its 
associated analytic signal when the used wavelet is a complex analytic func- 
tion (i.e. whose Fourier transform is zero for negative frequencies). Indeed, 
the single modulated components present in the signal may be extracted 
(provided they are asymptotic and sufficiently separated in frequency) from 
the restrictions of the wavelet transform to the curves where the wavelet 
modulus is maximum, which are called the ridges of the transform. More 
explicitly, if the signal is of the form 



X{t) = ^^=l^n{t) = K=lMt) COs[(^„(t)] = K=lMt) COsK(t) ■ t + y^on] 



1 d(p^{t) 
2n dt 
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and a Morlet wavelet ijjit) = e^'^ot^-— jg y^^Q({^ then the restrictions of the 
wavelet transform to the N ridge curves a = ttmit) (where the wavelet scale 
a is proportional to the inverse of the frequency) may be expressed as 

W^niarnit),t) = Corr{t)Z-,n + r{t), 

where the term Corr{t) is completely defined by the used wavelet and by the 
ridge values and the residual r(t) is negligible if the component is asymptotic. 
It is then possible to approximately obtain the analytic signals associated to 
the single components and to estimate the relevant frequency and amplitude 
modulation laws. However, some limitations exist in the frequency modula- 
tion laws that may be detected using this procedure [22] • In particular, we 
note that the signals that do not fit satisfactorily the asymptotic condition 
make the wavelet ridge analysis much less accurate than the usual Hilbert 
demodulation technique, see [THj. Consequently, a new procedure was de- 
vised, which in a sense exploits the qualities of both techniques and which 
will now be briefiy described. The signal is first transformed using a Morlet 
wavelet and the ridges are approximately derived from the modulus maxima 
curves; in this procedure, a high value of the central frequency of the wavelet 
uo is used, in order to assure a high frequency resolution and the reduc- 
tion of the interference effects between adjacent components. The wavelet 
maps are then filtered, neglecting the coefficients outside a band around the 
dominating ridge, and an inverse transform is applied. The extracted sig- 
nal may then be subtracted from the original one and the procedure may 
be repeated until all the detectable components are extracted. The Hilbert 
transform technique is then applied to each component, in order to obtain 
the required modulation laws. In the latter step a further filtering may be 
used, considering only the time intervals in which the modulus of the asso- 
ciated analytic signal is higher than a given threshold (in terms of a given 
percent of its mean value). This avoids spurious large fiuctuations of the 
instantaneous frequency and the final statistical analysis may be restricted 
to those intervals of time in which this frequency is physically meaningful, in 
the sense that it may be confidently associated with fiuctuations produced 
by a fiuid dynamical mechanism, like the passage or oscillation of vorticity 
structures. Obviously, in order to characterize their level of chaoticity, the 
single components obtained through the wavelet filtering procedure may also 
be analyzed with different techniques, as for instance those described in the 
following sections. 
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3 Diffusion Entropy Method 



This technique of analysis of time series was introduced originally in the 
article of Ref.^j, whose authors coined in fact the term Diffusion Entropy 
(DE). It became clear immediately afterwards that the DE method is 
an efficient method of establishing the scaling value if the scaling condition 
applies, regardless of the scaling origin, and even when the scaling condition, 
according to the Generalized Central Limit Theorem (GCLT) yields a 
divergent variance. Given a discrete signal we interpret it as a 

set of diffusion generating fluctuations. The collection of these fluctations 
yield a single diffusion trajectory. We can we create a set of many different 
diffusion trajectories by means of moving windows of size t with 1 < t < A^. 
We generate N — t + 1 trajectories considering the sums: 

i=j+t 
i=j 

Each diffusive trajectory can be thought of as the final position of a walker 
which jumped for t time steps. Let p{y, t) be the probability to be at position 
y after t time steps. If the scaling condition applies to the asymptotic time 
limit, p{y, t) is expected to fulfill the following condition: 

Ml/,t) = ^i^(|). (3) 

Some of the processes of anomalous diffusion fit the prediction of the GCLT, 
thereby assigning to the function F{y) the structure of a Levy distribution 
with a diverging second moment. This is a problem for the techniques of 
analysis based on the observation of the second moment. It is not a prob- 
lem for the DE, which evaluates the scaling parameter 6 from the Shannon 
entropy of the probability p{y, t): 

/ + 00 
p{y,t)\n{p{y,t))dz. (4) 
-co 

Indeed, it is straightforward to prove that when the scaling condition of 
Eq.(jni) holds, Eq.(jH) yields S(t) = A + 6ln t, where A is the Shannon entropy 
of the generating function F. Thus, the scaling parameter 6 is easily evalu- 
ated by plotting the Shannon entropy in a diagram with linear ordinate and 
a logarithmic abscissa. The Shannon entropy, as a function of time, becomes 
a straight line whose slope is the scaling parameter 6. 

It is important to point out that in the case where the time series under 
study is a periodic process, the diffusion entropy, starting from a small initial 
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value, reaches a maximum, then it regresses to the initial value, at a time cor- 
responding to the period of the process under study |7] . Then, moving from 
this value, the diffusion entropy increases again up to a maximum which has 
the same value as the earlier maximum, and so on. In the case of the logistic 
map at the chaos threshold the diffusion entropy exhibits a behavior very 
similar to that of a periodic process. The maxima still lie on an horizontal 
line, but the regressions to the minima are not complete as in the periodic 
case and occurr in a disordered way [Zj. 

In the case of periodic processes the initial entropy increase is due to the 
fact that the moving window technique, adopted to create a large number 
of different diffusion trajectories, is equivalent to setting uncertainty on the 
initial condition. 

4 Computable Information Content Method 

The second method of analysis used in this paper aims at establishing a direct 
contact with Algorithmic Information Content (also known as Kolmogorov 
complexity). In this second method, the basic notion is the notion of infor- 
mation. Given a finite string s (namely a finite sequence of symbols taken in 
a given alphabet), the intuitive meaning of quantity of information Iaic{s) 
contained in s is the length of the shortest binary message from which we 
can reconstruct s (e.g. if we are working on a computer, that binary message 
is a program p that outputs s) . This concept is expressed by the notion of 
Algorithmic Information Content (AIC). We limit ourselves to illustrating 
the intuitive definition of the concept of AIC; for further details see and 
[T7j and related references. 

By definition, the shortest program p which outputs the string s is a sort 
of optimal encoding of s: the information that is necessary to reconstruct 
the string is contained in the program. Unfortunately, this coding procedure 
cannot be performed on a generic string by any algorithm: the Algorithmic 
Information Content is a quantity which is not computable by any algorithm 
(see Chaitin theorem in |25j). 

Another measure of the information content of a finite string can also be 
defined by a loss-less data compression algorithm Z satisfying some suitable 
optimality properties which we shall not specify here. Details are discussed 
in ^Tj. We can define the information content of the string s as the binary 
length of the compressed string Z{s), namely. 



The advantage of using a compression algorithm lies in the fact that, this 




(5) 
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way, the information content Iz (s) turns out to be a computable function. 
For this reason we shall call it Computable Information Content {CIC). 

The notion of information is strongly related to chaos, unpredictability 
and instability of the behavior of dynamical systems. The KS entropy can 
be interpreted as the average measure of information that is necessary to 
describe a step of the evolution of a dynamical system. 

We have seen that the information content of a string can be defined 
either with probabilistic methods (following the Shannon theory) or using 
the AIC or the CIC. Similarly, also the KS entropy of a dynamical system 
can be defined in different ways. The probabilistic method is the usual one, 
the AIC method has been introduced by Brudno [2^]; the CIC method has 
been introduced in |2ZI and [201 ■ principle, it is possible to define the 

entropy of a single trajectory of a dynamical system. There are different 
ways to do this (see [SE], [2H], |Hj, [Oj, [HI!)- In this paper, we make use of a 
method which can be implemented in numerical simulations. Now we shall 
describe it briefly. 

Through the usual procedure of symbolic dynamics, given a discrete in- 
finite trajectory x = |xj = T*(xo)}^^q drawn from the dynamical system 
{X,fi,T), we consider a finite partition V = {Ri, . . . ,Ri) of the dynamical 
system. In a standard way, we associate a string $p (x) to the trajectory x. 
It is that $-p (x) = {sq, Si, . . . , Sk, ■ ■ ■) if and only if 

W k >0 Xfc G Rs^ , where G {1, . . . , /}. 

We can define the Information content I{x, V, n) of the trajectory x with 
respect to the partition V in the following way 

I{x,V,n) :=/(<l>p (xD , 

where $p (x)"' is the string made of the first n digits of the symbolic trajectory 
$p (x). The Information Content I can be measured either via AIC or via 
CIC, so we have just defined Iaic and Iz (respectively). 

Let us assume that the compression algorithm Z is optimal in the sense 
of Ref. [Ej. We have the following results (fZj). 

Theorem 1. If Z is an optimal coding, {X,fi,T) is an ergodic dynamical 
system andV is a measurable partition of X , then for /i- almost all trajectories 
X drawn from the dynamical systems it holds: 

Iz{x, V, n) = Iaic{x, V, n) = nh^{T, V) + o{n) , 

where h^{T,V) is the Kolmogorov entropy of {X,fi,T) with respect to the 
measurable partition V. 
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Theorem 2. Given the dynamical system {X,T,n), if the measure fi on 
X is T-invariant, then, if Z is an optimal compression algorithm, for any 
measurable partition V it holds 

h^{T,V) = [ limsup — \ djj, = f limsup — \ _ 

J X n^+oo n J X n^+oo ^ 

Let us set lim sup^^^j.^^ H^^^^^ri) ^j^^ complexity K{x,V) of the trajec- 
tory X with respect to the partition V. Theorem |21 shows that if a system 
has an invariant measure, its entropy with respect to a given partition can 
be found by averaging the complexity of its orbits over the invariant mea- 
sure. Then, the entropy may be alternatively defined as the average orbit 
complexity. However, if we fix a single point, its orbit complexity is not yet 
well defined because it depends on the choice of a partition. We chose to 
get rid of this dependence by considering only a particular class of partitions 
and define the orbit complexity of a point as the supremum of the orbit 
complexity over that class. 

Let Pi be a family of measurable partitions such that lim diam{j3i) = 0. 

If we consider the quantity limsup Kz{x,f3i) we have that the following 

i — >oo 

lemma holds (for the proof, see [T7j). 

Lemma 1. // (X, /i, T) is compact and ergodic, Z is optimal, then for /x- 
almost all points xq & X , for the trajectory x of starting point xq it holds 

limsup Kz{x,l3i) =limsup KAic{x,(3i) = h^{T) . 

i— >oo i— >oo 

Therefore, this lemma, under suitable conditions on the dynamical sys- 
tem, on the partitions and on the compression algorithm, permits to define 
almost everywhere the complexity of the orbit x of the starting point xq as 
limsup Kz{x,Pi) =limsup Kaic{x, Pi). This quantities asymptotically ap- 

i— »oo j— »oo 

proach the KS entropy of the system. Indeed, the asymptotic behavior of 
I{x,V,n) gives an invariant of the dynamics which is finer than the KS en- 
tropy and is particularly relevant when the KS entropy is null. 

It is well known that the KS entropy is related to the instability of the 
orbits. The exact relations between the KS entropy and the instability of the 
system are given by the Pesin theorem. We shall recall this theorem in the 
one-dimensional case. Suppose that the average rate of separation of nearby 
starting orbits is exponential, namely, 

Ax{n) ~ Ax(0)2^" for Ax(0) < 1, 
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where Ax{n) denotes the distance at time n of two points initially at distance 
Aa;(0). The number A is called Lyapunov exponent; if A > the system 
is unstable and A can be considered a measure of its instability (or initial 
data sensitivity). The Pesin theorem implies that, under some regularity 
assumptions, A equals the KS entropy. 

In weakly chaotic systems, the amount of information necessary to de- 
scribe n steps of a trajectory is less than linear in n and the rate of separation 
of nearby starting orbits is less than exponential, then the KS entropy is not 
sensitive enough to distinguish the various examples of weakly chaotic dy- 
namics since in every case the KS entropy is zero. Nevertheless, using the 
ideas we illustrated above, the relation between initial data sensitivity and 
information content of the orbits can be extended to these cases. 

An example of such a generalization is given in Ref . • In the following, 
for the sake of brevity we only briefly sketch the ideas underlying the results 
of Ref. fO] , which are very deep and detailed. Let us consider a dynamical 
system ([0,1], T) where the transition map T satisfies some constructivity 
properties (a constructive map is a map that can be defined using a finite 
amount of information), and the function I{x, V, n) is defined using the AIC 
in a slightly different way than before (use open coverings instead of parti- 
tions). Provided the speed of separation of nearby starting orbits goes like 
Ax{n) ~ Ax{0)f{x,n), it has been proved in Ref. [lSi\ that under suitable 
assumptions on /, we have for almost all the points a; e [0, 1] 

/(x,P,n) ~ log(/(x,n)). (6) 

In the weakly chaotic case, the speed separation of nearby starting or- 
bits is less than exponential. In particular, if we have power law sensitivity 
Ax{n) ^ Ax(0)nP, the information content of the trajectory is 

I{x,V,n) ~plog(n) . (7) 

If we have a stretched exponential sensitivity Ax{n) ~ Ax(0)2^"*', p < 1 
(e.g. in the Manneville family with z > 2) the information content of the 
orbits will increase with the power law: 

I{x,V,n)^nP. (8) 

Since we have shown that the analysis of I{x, V, n) gives useful informa- 
tion on the underlying dynamics and since J(x, P, n) can be defined through 
the CIC methodology, it turns out that it can be used to analyze experimen- 
tal data using a compression algorithm which is efficient enough and which is 
fast enough to analyze long strings of data. For the numerical experiments. 
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we have used a compression algorithm called CASToRe, illustrated in an 
earlier publication [20]. The name of the algorithm is an acronym meaning 
Compression Algorithm, Sensitive to regularity; the heuristic motivations for 
this name are explained in Therefore, in the following the algorithm Z 
will be the algorithm CASToRe. 

5 Results 

As a significant earlier work on real fluid dynamical processes, we have in 
mind the pioneering work of Solomon, Weeks and Swinney [22^ , who studied 
chaotic transport in a laminar fluid flow in a rotating annulus by tracking 
large numbers of tracer particles for long times. The remarkable result of 
that paper is the detection of Levy statistics. This means that these authors 
discovered strange kinetics. Our approach is slightly different. We do not 
adopt the Lagrangian approach to chaotic transport of Ref. [201, but rather 
the Euler perspective: we measure the velocity fluctuations in a suitable 
region of the flow. 

In a jet flow, the shear layer between the so-called potential core and 
the ambient fluid is subjected to an instability, which produces a roll-up of 
the shear layer into vortical structures; moving downstream, these structures 
grow and normally undergo a process of pairing, through which two subse- 
quent structures merge to form a larger one. Depending on the geometrical 
and fluid dynamical conditions, multiple pairings may occur before azimuthal 
instabilities lead to a complete mixing of the flow at the end of the jet core. 
The instability leading to roll-up is characterized by a dominating frequency, 
while the pairing process leads to a doubling of the period of passage of the 
structures, and then to the appearance of its subharmonics. In the case of 
a coaxial jet (Figure HJ, the situation is complicated by the presence of two 
cores and shear layers, which give rise to different families of vortices whose 
degree of interaction is a function of the ratios yf- and jf-. In the present 
case = 2,^ = 1.5), the stronger outer shear layer dominates the initial 
development of the flow and its dynamics is not far from that of a single jet 
PT] . In order to describe the performance of the different signal processing 
techniques in providing information on this dynamics, in the following we 
shall analyze the axial velocity signal u obtained from the measurements in 
the position = 1, = 0.96, i.e. well inside the outer shear layer. 

As shown in Figure El the Fourier spectrum of this signal has two peaks, 
one at the fundamental frequency near 190 Hz (corresponding to a Strouhal 
number, based on Uo and on the outer shear layer momentum thickness at the 
jet outlet, St = f9/Uo ~ 0.0125), and a larger one at the subharmonic. This 
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indicates that at this position the pairing process is already in an advanced 
stage, although not yet completed, as may be derived from the analysis of the 
radial velocity spectrum (not shown here for the sake of brevity), in which 
the peak at the fundamental frequency is still the dominating one. 

The two components of the u signal were extracted using the wavelet 
filtering procedure described above and we reconstructed the two signals Ua 
and Uh, associated to the harmonic 190 Hz) and to the subharmonic {ub, 
95 Hz) frequencies, respectively. The two reconstructed signals Ua and Ub 
were then analyzed using the different techniques described in the previous 
sections, starting from the Hilbert demodulation procedure. 

Figure El shows a portion of the time variation of the instantaneous fre- 
quencies of the two components derived from their associated analytic signals. 
As already pointed out, to avoid unphysical variations of this frequency due 
to the degeneracy of the phase definition, the data are filtered considering 
only the time intervals in which the modulus of the analytic signal is above 
a given threshold (50% of the modulus mean value in Figure Ej) and are then 
statistically analyzed. By this procedure, the mean values of the two frequen- 
cies were found to be 188.80 Hz and 94.47 Hz, with standard deviations of 
4.8 Hz and 0.93 Hz, respectively. This technique also allows us to obtain the 
time variation of the ratio between the two instantaneous frequencies, whose 
mean value was found to be 2.002, with a standard deviation of 0.05, thus 
confirming that Ub is the subharmonic of Ua- As can be seen from Figure El 
the instantaneous frequencies vary in time with an irregular trend, and the 
same may be shown to be true for the modulus of the analytic signal, which 
corresponds to the amplitude of the original velocity component fluctuation. 
This irregular behavior is most probably due to the fact that the shape, size, 
radial position and translation velocity of the vortical structures are not ex- 
actly constant in time. The nature of this irregularity will now be analyzed 
through the DE and CIC methods. 

First, we have studied the two reconstructed signals Ua and Ub using the 
DE method. The results are illustrated in Figure IH and, to make our dis- 
cussion more transparent, are compared to the result provided by the same 
method in the case of a periodic signal. The results of Fig. ^ deserve some 
detailed comments. First of all, we note that the two reconstructed signals 
yield a similar behavior, with an initial transient followed by an oscillatory 
behavior with the DE maxima lying on an horizontal line. This is property 
shared by the periodic signal. It is evident, however, that the similarity with 
the periodic signal is not complete. Indeed, in the case of the periodic signal 
the regressions of the DE to the initial condition are complete, while in the 
case of the two reconstructed signals they are not; in fact, the regression 
intensities tend to decrease with time, suggesting that at infinite time the 
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Figure 2: Fourier power spectrum of the axial velocity u at position x/Di 
and z/Di = 0.96. 



diffusion entropy might become rigorously constant. This perfect locahzation 
is similar to the slow occurrence of a collapse noticed year ago in the case 
of a problem of quantum chaos [30^. In that case the dynamics under study 
was that of a 1/2 dipole precessing about a magnetic field while undergo- 
ing the perturbation stemming from the harmonic motion of an oscillator. 
The entropy recursions to the initial condition are complete if the Larmor 
frequency is very small compared to oscillator frequency. In the case of re- 
sonance, however, chaos emerges and consequently incomplete regressions of 
decreasing intensity were recorded. All this suggests that the process un- 
der study is neither periodic, nor quasi-periodic. In fact, in both cases the 
entropy regressions to the initial condition would be complete |7| IHlj. 

In order to evaluate the degree of regularity of the two reconstructed sig- 
nals, we have estimated the dominant Lyapunov exponent of the two time 
series [32|, which turned out to be zero for both sequences Ua and Ub- There- 
fore, thanks to the Pesin Theorem, we can say that the vortex dynamics in 
both cases of roll-up and pairing has null KS entropy. Together with the fact 
that the processes are not ordered, we conclude that the vortex roll-up and 
pairing are weakly chaotic phenomena. 

As a consequence, we classified the kind of complexity of the vortex dy- 
namics by investigating the behavior of the information content of the recon- 
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Figure 3: Instantaneous frequencies of the components Ua andub , as obtained 
using the wavelet-Hilbert procedure. 



structed signals. We have translated the time series into symbolic sequences 
by considering the interval Ja = [mm{ua},Taax{ua}] (similarly, we defined 
Jb) and dividing it through several uniform partitions Vk with k subinter- 
vals where k G {2,4,8,16,32,64,100}. As a result, the different associated 
symbolic sequences gave the same power-law behavior for the information 
content, which has been calculated by means of the algorithm CASToRe: 



The exponent a is 0.80 for Ua and 0.83 for Ub, in both cases it does not 
depend on the particular partition used. The constant C depends on the 
partition Vk and increases with the parameter k (see figure E]). As illustrated 
in section 4, the fact that the amount of information to describe n steps of a 
trajectory is less than linear in n is in complete agreement with the null KS 
entropy. 

At this stage, we have to address the challenging task of building up a 
model reproducing the main properties of vortex dynamics found so far. This 
model must explain the degree of complexity as it is indicated by the parame- 
ter a ranging from 0.80 to 0.83, and at the same time it has to be compatible 
with another crucial property such as the two-peaks Fourier spectrum of Fig. 
2. We propose a theoretical model, resting on the basic idea that the exper- 
imental signals are driven by a hidden intermittent process inspired to the 
Manneville map. We shall use this model to reproduce both the signal Ua 



Izin) ~ Cn' 



a 



(9) 



15 




Figure 4: results of the DE analysis on: Ub (top graph, dashed line), Ua (mid 
graph, solid line) and a purely periodic signal (bottom graph, dotted line). 
The scale is log-linear. 



and the signal Uh- However, we shall give details only on the choices made 
to reconstruct Uh- The prescriptions adopted to reconstruct the signal Ua are 
similar, and, for the sake of simplicity, will not be illustrated. First of all, we 
define the following waiting time distribution: 

= - . (10) 

Note that this form is chosen as the simplest way to generate an inverse power 
law at long times, without conflicting with the normalization condition. In 
the case /x > 2, the first moment of this distribution is T/(/i — 2), thereby 
showing that the parameter T keeps under control the mean value of this time 
distribution. We also define the set of frequency values A(, = {ui, ....,u!e} 
where the values tu^'s are determined on the basis of the Hilbert transform 
approach used to derive the results of Fig. HI More precisely, these values 
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Figure 5: The information content of the signal Ui, as a function of the time 
length t, in a log-log representation. We show the results corresponding to 
different partitions , al of them sharing the same power law exponent 0.83. 
From top to bottom: partition Vioo, V32, T^S; V2-The bottom solid line is the 
analytical prescription for the power-law exponent 0.83. 



can be considered to be a coarse-graining representation of the condition 
illustrated in Fig. IHl In the specific case of the signal Ub, we have fixed i = 10 
and the frequencies cjj's have been determined by a uniform partition of the 
frequency centered at (94.47 Hz) with a width that is twice the standard 
deviation (0.93 Hz). 

We now define the artificial signal as follows. We select a time ri of 
the distribution of Eq. pO|l . we call laminar phase this time interval and we 
randomly select one of the frequencies of the set Ai, say lj{ti). We assume 
that the signal, a sinuisoidal function of time, keeps this frequency throughout 
the whole laminar region. At the end of this laminar region, we select another 
number T2 from the distribution iP{t) and another frequency uj{t2) from the 
set Ai, and so on. At any random drawing of the pair {xj, ci;(rj)}, we must 
also select a phase (pj. This latter choice is not arbitrary, but it is done in 
such a way as to ensure that the resulting signal is continuous. In conclusion, 
we get the following signal 
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Figure 6: result of the DE analysis on: Ub (top graph, solid line) and the sim- 
ulated signals cr^(t) generated with different values of the control parameter. 
From top to bottom: fi = 3.5, /i = 2.45, fi = 1.90. 



where the subscript /i indicates the dependence of this artificial signal on 
the random distribution of Eq. (fTUj) . with a fixed value of /x. This special 
way of sewing a laminar region to the next has been dictated by the need of 
reproducing the experimental spectrum of Fig. 2. Our numerical calculations 
allowed us to assess that a different choice, based on the abrupt frequency 
jump , with a fixed phase, and consequently with an abrupt jump of the signal 
(Tfj,{t), yields a continuous Fourier spectrum, thereby dramatically departing 
from the behavior illustrated by Fig. 2. 

The inverse power law distribution given by ()1()|1 is an idealization of 
the Manneville map mentioned in the introduction. Nevertheless, the infor- 
mation content growth of the instantaneous frequency ujfj,{t) is the same as 
that of a trajectory drawn from a Manneville map with driving parameter 
z = -j^. Therefore, the three types of statistics illustrated in the Introduc- 
tion become: /i > 3, ordinary statistics; 3 > /i > 2 strange kinetics in the 
sense of Zaslavsky P] and fi < 2 sporadic dynamics in the sense of Gas- 
pard and Wang J5]. The sewing process adopted to reproduce the Fourier 
spectrum of the real signal has the effect of bidding this complexity so as to 
make the DE analysis insensitive to the complexity of this hideen process. 
This is clearly illustrated in Fig. 6, where we see that the artificial sequences 
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with three different degrees of hidden complexity, ranging from weak chaos 
(/i = 1.90) to ordinary statistical mechanics (/i = 3.5) through a condition 
of strange kinetics (/i = 2.45), result in the same time behavior of S{t). We 
note that, in spite of a lack of sensitivity to the complexity of the hidden 
driving process. Fig. IHl proves that there is a good agreement between the 
DE of the reconstructed signal Ub and the DE of three simulated sequences 
generated with different driving parameter /i. Therefore, Eq.illlj). with the 
random choice of waiting times prescribed by Eq. (|TO|) . is a plausible model 
of the turbulent process under study. 

In conclusion, while the complexity of the reconstructed signals is mea- 
sured with the compression algorithm and turns out to be given by Eq.Q 
with a = 0.80 in the case of and a = 0.83 in the case of Ub (see Fig. the 
complexity of the driving process uj{t) remains hidden. At the moment we do 
not have any analytical expression to relate the information content index 
a to the driving parameter /i of Eq. (fTT| . Thus, it is impossible to deter- 
mine theoretically the degree of complexity of this driving process. We must 
rest on a numerical treatment, based on applying the compression algorithm 
to artificial sequences with different values of fi. For this purpose, we have 
turned the artificial time series generated by Eq. ljllj) into several symbolic 
sequences by means of a uniform partition V2- For all the values of fi adopted 
in the numerical simulations, the resulting information content of associated 
symbolic series turned out to be Eq.©, with a < 1, which is equivalent 
to the sporadic randomness of Gaspard and Wang |H]. The dependence of 
a on /i is pictured in Figure and turned out to be essentially monotonic. 
The numerical function a(/i) crosses the horizontal strip 0.80 < a < 0.83 
at /i = 2 ± 0.1. Thus the complexity of the driving stochastic process is 
measured, with a good approximation, by equation (fTT|) with fj, = 2. 

6 Concluding remarks 

The interest of this paper rests on the discovery of the hidden complexity 
of the signals Ua and Ub, extracted from the vortex dynamics by means of 
the wavelet analysis. The main properties of these signals are satisfactorily 
reproduced by the artificial signal of Eq. (fTTjl . obtained according to the 
procedure described in Sectional This important aspect is illustrated by Fig. 
|71 For this reason we devote these concluding remarks to further comments 
on this figure and we stress that from it two main results emerge. 

The first result is that the sewing process adopted to generate the artificial 
sequence yields a sort of time dilution of randomness. The driving signal for 
H > 2 is characterized a non- vanishing Lyapunov coefficient and consequently 
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Figure 7: The index a denoting the complexity of the signal cr^{t) of Eq. 
(10) as a function of the parameter fi, denoting the complexity of the hidden 
driving process. 



by an Information Content growing linearly in time. The time increase of 
the Information Content becomes sub-linear only in the region /i < 2. We 
see from Fig. [7|that the Information Content of the signal cr^{t) increases 
sub-linearly in time throughout the whole range of yU-values considered, and 
remains probably sub-linear well beyond /i = 3, which is the border between 
Levy s and ordinary Gaussian diffusion. According to the definition of weak 
chaos adopted in this paper, we conclude that the signal is an expression 
of weak chaos regardless of the degree of hidden complexity. We do not have 
available at the moment a theory expressing a as a function of /i. However, 
the numerical results of Fig. [3 indicate that a{n) is a monotonic function 
of /i with a < 1. The second result has to do with the detection of the 
degree of hidden complexity, namely, the complexity of the driving signal 
uj^{t). The uncertainty on the complexity of the signal u is measured by 
the width of the strip between the straight line a = 0.83 and the straight 
line a = 0.80. Thus, the uncertainty on the complexity of the hidden signal 
can be determined by observing the crossing between the experimental curve 
a{fi) and the strip. We see from Fig. El that the crossing occurs in a sharp 
interval of /i values around /i = 2. More precisely, the complexity of the 
driving signal is measured by 1.9 < /i < 2.1. This makes us conclude that 
the complexity of the driving signal lies at the border between weak chaos 
(yU < 2) and strange kinetics {2 < n < 3. There are other complex processes 
lying very close to this border (see, for instance, the solar flares of Ref. (SB] 
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and the earthquakes in Cahfornia It would be interesting to assess why 

the border between the stationary and the non-stationary states is the basin 
of attraction of some complex systems. 
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